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Abstract. We present a technique novel in numerical methods. It compiles 
the domain of the numerical methods as a discretized volume. Congruent 
elements are glued together to compile the domain over which the solution 
of a boundary value problem of a linear operator is sought. We associate 
a group and a graph to that volume. When the group is symmetry of the 
boundary value problem under investigation, one can specify the structure of 
the solution, and find out if there are equispectral volumes of a given type. 
We show that similarity of the so called auxiliary matrices is sufficient and 
necessary for two discretized volumes to be equispectral. A simple example 
demonstrates the feasibility of the suggested method. 



1. Problem description 

In both science and engineering, we solve boundary values in a volume composed 
of large number of meshes. This is the case in nuclear engineering [T], in fluid 
dynamics [5], [3] and electromagnetic fields [1]. We address the question: under 
what conditions are the solutions for two meshes identical, or, are the solutions 
transformable into each other by a simple rule? How to find the transformation 
rule? How can we find equivalent meshes? 

In the design and safety analysis of large industrial devices, calculational models 
are tested against experiments carried out on a small scale mock-up. This is the 
case with nuclear power plants [5], aeroplanes [6], and ships [7]. We would need 
a transplantation of the measured values to the geometry of the real scale device. 
Is there any hope of doing that exactly or have we to put up with approximate 
methods [5]? 

In the sequel, we assume operator A to be a linear operator defined in a finite 
domain V in and to commute with the symmetry group of the plane . As to 
its physical meaning, the authors had in their minds the Laplace operator occur- 
ring in several physical problems from electricity to quantum physics, the diffusion 
and/or transport operators as used in reactor physics with homogeneous material 
distributions. The boundary value problem considered in the sequel is of Cauchy 
type (solution is zero along the boundary of V) , but the authors are convinced that 
generalization to Neumann or third type boundary value problem is straightfor- 
ward. 

We suggest using discretized volumes. The domain V, over which the solution of 
the boundary value problem is sought is constructed by the following procedure. We 
choose an appropriate simply connected tile t, which in our cases is an n-gon and 
glue copies of the tile by their corresponding sides. Although in principle t is almost 
arbitrary [9], we confine the discourse to triangular shaped tiles. Since triangulation 
is a well known and widely used technique even in theoretical problems [10], this 
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is not considered as a limitation. The discretized volumes considered by us are 
always finite. A concise description of the structure of a discretized volume V is 
a finitely presented group G . Although in practical problems the applications 
of computational group theory are rather limited, the authors strongly hope for a 
steady development in both computational tools (software) and means (hardware). 
So, a discretized volume is described by the tile t and the group G. A further asset, 
a graph F is also defined. If copies ti and tj of tile t are interconnected by an edge 
a of t, graph vertices i and j are also interconnected by an edge s^. Analyzing the 
group G, the graph F, one can easily reveal basic properties of V. The main results 
of the present work are: 

(1) We provide an algebraic description of the discretized volume. By analyzing 
the group and graph associated with a given discretized volume, one can 
answer a number of questions. 

(2) Using that algebraic description, we can formulate conditions for equispec- 
tral volumes to exist. 

(3) We formulate a formal solution to the eigenvalue problem to be specified 
later. 

(4) We give conditions for two discretized volumes to be equipectral and we 
show that the eigenfunctions of equispectral discretized volumes are trans- 
formed into each other by a linear map. 

The structure of the manuscript is as follows. We define the discretized volume 
(DiV) in Section 2, along with the associated group theoretic assets, group and 
graph. In Section 3, we present the solution space Ly, the space of square integrable 
functions. By the discretized volume being glued copies of a tile t, the solution 
is decomposed into functions defined along t. This allows for the dot product 
to be applied to functions over different discretized volumes. Section 4 discusses 
the structure of the solution for the boundary value problem under consideration. 
Section 5 treats a well known example by the tools presented in parts preceding 
Section 4 of the manuscript. 

2. Algebraic Description of Discretized Volumes (DiV) 

In the present section, we study special planar domains over which we are going 
to solve eigenvalue problems. Throughout the present work we investigate cases 
where i C is an acute scalene triangle (see Figs. 1,2). 

Definition 2.1. A discretized volume V is composed of < oo (finite number) 
copies of tile t so that we glue copies of tile t to each other along a corresponding 
side. Those copies of t, which share a side are called adjacent. The shared side is 
called internal. When A^ > 1, every copy of t has at least one internal side, see 
Figs. 1, 2. 

Definition 2.2. We say that the discretized volumes Vi, Vi are equivalent if there is 
an isometry of the Euclidean plane which maps V\ into V2 . Equivalent discretized 
volumes are denoted as V\ ~ V2. 

We only remark here, that definition 12.21 does not distinguish the "warped pro- 
pellers" [T7j because they are obtained from each other by interchanging the dotted 
and dashed sides of t and that operation leaves t invariant. Our method for con- 
structing DiVs has immanent limitations. Depending on the tile t, we may tile 
out the entire plane, or, after a given number of gluing, we have to stop because 
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the next glued copy would intersect with an already existing copy. To clear that 
problem, we investigate the transformation rules of tile t under gluing. 

Definition 2.3. We label the corresponding side of the congruent copies of tile t 
in DiV V by q;,/3,7. Let be a local reflection on the side v oi t, i' ^ {a,/?, 7}. 
The image of t under is denoted by — ■ t. 

By gluing, we get a new copy, t,^, and gluing is applicable to the new copy as well, 
so we can define an operation among the Ti, transformations and that operation is 
the consecutive application of gluing. 

Definition 2.4. Let = T^^T^^ mean the following transformation: apply first 
T^^ and then T^^ to the result: T^^T^J = T,,^ (T^J i, see Figs. 1, 2. 

We are going to use the term set of gluing to the set of transformations defined 
above. On the set of transformations we defined multiplication, and there is a unit 
element Te which leaves t invariant. Such an element is, for example, the repeated 
application of T^. Trivially there is also an inverse, because if is obtained by 
gluing from t then also t is obtained from ti, by gluing. 

Definition 2.5. The set of transformations {T^^^T^^ . . . T^^} is endowed with the 
multiplication, see definition 12. 4i there is an identity element Tg = T^T^, and 
every element has an inverse, therefore the transformations form a group Gj . 

We are going to refer to the group of transformations as Gf, where t refers to 
the tile, the corner stone of DiVs. Note that we allowed for a repeated application 
of a given transformation and then we get back the original tile. But there are tiles 
which after a sequence of transformations only partly coincide with a former copy 
of tile t. We exclude such a situation from our investigation. 

Definition 2.6. We call T eGf realizable if for any factorization of T = T^T^ 
the images • t and T^, ■ t are either disjoint or coincide for any T e G(. 

The set of realizable transformations depends solely on the tile t. From now on, 
we deal solely with realizable T transformations. 

Now we pass on to investigate the discretized volume V. We wish to associate 
algebraic descriptions with V, we define a graph and a group. 

Definition 2.7. A graph Fy is assigned to V in the following way. We label the 
copies of t in V. If the copies labeled as ni and n2 are adjacent, and they share a 
side of type a, then the vertices ni and n2 of graph Fy are connected by an edge 
of type a. 

Definition 2.8. We associate a permutation group G to in the following way. 
When in V side a of type a connects the copies U^^, ti^^^; tj^^, tja.2i ■ ■ ■ then, we 
form the permutation a = (iqi, 1^2), (jai, ja2)- • ■• We repeat that procedure for 
sides a, /3, 7 to get generators a, 6, and c, and group G is generated by a, 6, and c. 

The next step i to define the group action of G on V , and t, see Fig. 1, and Table 
1. There is a natural map of an element g gG to an element GGj. Either group 
is finitely presented, the number of generators is the same for both groups. Let 
G=< a, 6, c|a2 = 6^ = = e > and G* =< T,,, T^, T^IT^ = T^^ ^ rp^2 ^ ^ 
Then the map a — f{a),P = f{b), 7 = /(c) maps G— > Gt and is an injection. Such 
a permutation representation is a so called faithful representation. We apply the 
notation v — f{g), g £ {a, 5, c}. 
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Table 1. Group action in Fig. 1. 



generator/copy 


1 


2 


3 


4 


a 


2 


1 


3 


4 


b 


1 


3 


2 


4 


c 


1 


2 


4 


3 



Definition 2.9. The action ((7 • i) of g G G on a tile t is defined as T^, ■ t where 

Definition 2.10. The orbit of t under the group G is the set Tg ■ t for aU g G G. 

Definition 2.11. The action g ■ V oi g G G on 1/ is defined as follows. Let 
X G {a, 6, c} be a generator of G. The action of x on copy ti is x ■ ti = ti whenever 
side f~^(x) of copy ti is not internal in V . Otherwise x ■ ti = tj if copies tj and ti 
< N) share an internal side of type f^^{x). 

Remark 2.12. With a given tile t, N may be limited, see definition 12.61 

Definition 2.13. Adjacency matrix Ay of F is an x iV matrix, its aij element 
is aij = 1 if copies ti and tj are adjacent, otherwise a.y = 0. 

Definition 2.14. The auxiliary matrix X is X=D+ Ay, where £) is a diagonal 
matrix, its i*''cntry Da equals the number of internal sides of copy ti in V. 

Discretized volumes are shown in Fig. 2, with tile t, a regular triangle; with two, 
three, and four copies of t. The side types are solid (side a), dashed (side b), and 
dotted (side c) lines. The discretized volumes are described by t (identity element), 
{t^a -t}, {t,a ■ t, f3 ■ a ■ t}, {t, a ■ t, f3 ■ a ■ t,"f ■ P ■ a ■ t}. Let us number the copies of 
t as follows: t ^ 1, a ■ t ^ 2,/3 - a- t^S, '^■/3-a-t^4. The group action on the 
four elements is given in table [1] We have elaborated the basic means to be used 
to analyze the solution of an eigenvalue problem over a DiV. It is well known that 
it suffices to solve an eigenvalue problem over one element of the set of equivalent 
discretized volumes. This is because the equivalence provides a map between equiv- 
alent volumes and when that map is built from transformations commuting with 
the operator in the eigenvalue problem we immediately get a transformation of the 
solutions. We need further means to recognize if two complex DiVs are equivalent. 

Lemma 2.15. IfVi ~ V2 and the corresponding graphs are Ti and T2, then Fi ~ 

r2. 

Proof: If the stipulated conditions are met then there is an isomorphism between 
the copies and edges of the two discretized volumes this entails the statement. □ 

Lemma 2.16. Let the number of copies oft in discretized volume V be N. Then, 
in group G associated with V , there is a subgroup of index N . Consequently, the 
order of group G is a multiple of the number of copies N in V . 

Proof: There is a map between finite groups and Cayley graphs. A coset repre- 
sentation of G is isomorphic with graph F, in which there are A^ vertices. □ 

In accordance with Lemma [2.151 Gi ^ G2 suffices for Vi ~ V2. In accordance 
with Lemma [2.16| Gi ~ G2 suffices for Vi ~ V2. 

The next section deals with functions defined over a discretized volume. 



5 



3. Function space L\/ 

The specific structure of the DiV can be exploited in the analysis of eigenfunction 
defined on V. Let us consider a functioi0 f{x),x £ V. We trace back f{x) to A^- 
tuples (here N is the number of copies in V), which is a vector space. Thus we 
can speak of linear independence of two functions or the dimension of f{x). We set 
forth the following following notation. 

X = {xi,X2) denotes a point in the discretized volume V. In tile t, we use a 
local coordinate ^ = (^i,'C2)- Since V is composed of copies of t, and copy ti is 
obtained as ti = T; • t, and transformation t — > is the automorphism of the plane, 
in other words a member of the Euclidian group E(2), the following coordinate 
transformation, which acts on triples {xi, X2, 1) and is associated with 1^: 

(cos 9i sin 9i 
sin cos 9^ 
Qi bi 1 

That transformation maps a point x = {xi,X2) into 

(3.2) xgi — {xi cos9i + X2 sin9i + ai, — xi sin^^ + X2 cos9i + bi). 

The map has the property xgigj — {xgi)gj, gi, gj G E(2). Action of gt on a function 
f{x) is 

(3.3) TJix) - fixg,). 

Using the above definitions, we have a transformation ti ^ t as 

(3.4) e-2:.gri. 

For a given x € V also holds x G ti for some 1 < i < A^, i.e. point x belongs to one 
of the copies of tile t. This gives rise to a map f{x) to /(^) a value at point ^ of t: 

(3.5) iT,)-'f{x)^MO,xetf,^et. 

When V is composed of N copies of tile t, any f{x) is exhaustively described by the 
TV-tuple of functions /(O = (/i(0, /2(0, /jv(e)), where /,(0 = {T^r\nx),x G 
t^. 

Definition 3.1. The TV-tuple / = (/i(e), /2(e), • • • ,/iv(0), where /.(^ = iT,)-^f{x),x e 
ti is called the vector form of function f{x) defined over the discretized volume V. 

Definition 3.2 (Function space Ly). The function space hv contains square in- 
tegrable functions over discretized volume V . 

Definition 3.3 (Dot product and norm). Let f,h G hy- Then the dot product of 
/ and h is 

N 

(3.6) {f,h)v = Y. mum, 

where /i(^) = {Ti)^^ f{x),x G ti. The notation for hi is analogous. The norm of 
f{x) G Lv is 

N 

(3.7) \\f{x)f = Y. //'m- 



"'^We use the notation f{x) as a general function, not necessarily the same as before. 
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It is evident that the dot product ()3.6|) meets the foUowing general properties of 
the dot product (here /, /i, /2, h are functions, ai and 02 are numbers): 

• symmetry: 

(/, h)v^{hj)v; 

• Schwartz inequality: 

(/,Mi^<ll/lll|/^ll; 

• linearity 

(ai/i +a2/2,/i) = ai{fi,h)v + 02(72, /i). 

Definition 3.4 (Dimension of f{x)). dim{f{x)) is the number of linearly indepen- 
dent /i(^) functions in the vector form of f{x). dim{f{x)) is called the dimension 
of function f{x), defined over the discretized volume V. 

The dot product p.6p is actually formulated through the vector forms of the in- 
volved functions: 

N 

(3.8) (fix), h{x))v = itiOgm = E / M0h^{0d^■ 
Since 

(3.9) iMl+h)t - U+M^h)t 

for any N x N matrix M , we may use the usual rules of vector dot products: 

(1) imAmMUm,UhiO)t if f/ is a unitary matrix. 

(2) {M^l,h)t^{LMh)t 

(3) vectors / and h are called orthogonal if {f,h)t — 0, or, what is equivalent, 

if (/(x),l(x))y =0. 

Note that solely N, the number of copies in the volume of the integration, and the 
tile t are relevant in the dot product, therefore / and h may belong to different 
discretized volumes, provided each of them is composed of N copies of tile t. 

The next section deals with the structure of the solution to the eigenvalue prob- 
lem. 

3.1. Equispectral discretized volumes. The goal of our investigation is to find 
out if there are non-equivalent DiVs that allow for transforming the solutions into 
each other. We follow Baron Miinchausen's procedure[|: we assume that the solution 
is known along the internal boundaries of the DiV and give a formal solution in 
term of that. This leads us not only to the structure of the solution but also to a 
linear transformation mapping the solutions over two DiVs into each other. 
Let us investigate the following eigenvalue problem in a DiV Vi: 

(3.10) A^{x) = X^x), x G Vi. 

We assume ^{x) to be in C^(Vi), i.e. the function and its first derivative are 
continuous in V . Furthermore we assume A to commute with the automorphisms 
oiVi: 

(3.11) [A,0] = AO - OA = 0. 

Here O is an automorphism of T^, viz. a reflection, translation or rotation operator 
acting on functions defined over Vi . Eq. (|3.1ip is the usual definition of symmetries 



'Tale hero Baron Munchausen once fell into a swamp and drew out himself by his own forelock. 
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of operator A. The set of operations with which Vi is formed involves only symme- 
tries of operator A. Definitions 12.91 and 12.111 assure group G, defined there, to be 
isomorphic to a group of transformations commuting with A. It is well known that 
the solution of eigenvalue problem (|3.10p with a homogeneous boundary condition 
along the boundary dVi of Vi is easily determined from a solution of the same 
problem over another member of the class ~ Vi (see Definition 12. 2p . The question 
is, if we can find a volume V2 not equivalent to Vi such that all the eigenvalues of 
problem (|3.10p will remain the same as for Vi . 

Definition 3.5. If there exist volumes Vi and V2 each one composed of the same 
number of copies of the same tile t, so that Vi and V2 are not equivalent according 
to definition 12.21 furthermore all the eigenvalues of problem (|3.10p are the same in 
Vi and V2, we call Vi and V2 equispectral. 

When speaking about a fixed operator A, its spectrum is different on different 
DiVs. We express the formal solution in terms of assumed known solution along 
internal boundaries. Can those functions along an internal boundary be identically 
zero? Below we show that either the solution along each internal boundary is 
identically zero-that is the degenerate case-or, on every internal boundaries the 
solution differs from zero. 

Definition 3.6. The set S of eigenvalues A in Eq. ()3.10|) supplemented with ^(.to) = 
0, xq £ dV is called the spectrum of operator A on DiV Vi. 

Lemma 3.7. Let Svi be the spectrum of operator A on DiV Vi composed from 
tile t. Let St be the spectrum of operator A on t. Let Ao 6 5** H Svi and $0(2:) be 
the associated eigenfunction in Vi . Then $0 (x) is identically zero on every internal 
boundary of Vi . 

Proof: See [12]. □ 

The benefit from knowing equispectral volumes comes from the fact that it is 
rather tiresome to solve Eq. (|3.10p even for a simple volume. At the same time we 
know that equivalent volumes provide an easily feasible recipe for transplanting the 
solution from one member of the class to another. Knowledge of equispectral vol- 
umes would widen the range of transformations where instead of solving Eq. (|3.10p 
over a new equispectral volume, one would apply a relatively simple transformation 
to an already known solution. 

Actually, in a number of practical problems, one would be satisfied with the 
equivalence of one eigenvalue. In other words, with the equality of the respective 
eigenvalue and a transformation rule for the eigenfunctions. In a number of cases 
physical meaning is attributed to the so called fundamental mode eigenvalue. 

4. Solution of a Boundary Value Problem over DiVs 

4.1. The formal solution. Below we derive a formal solution [22] to problem 
(|3.10p with Dirichlet boundary condition (i.e. $ = on the boundary). The 
solution is given in terms of the Green's function Qt of tile t, that we obtain as the 
solution of the following boundary value problem: 

(4.1) (A-A)gt(e,xo) -o,eet, 



(4.2) 



Gt{^,xo) ^ 6{£, ~ xo), xo e dt. 
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Our goal is to express the solution of ()3.10p in terms of given values along the 
boundaries. In general, the boundary value uniquely determines the solution. In 
order to build up the solution in DiV V, we build up the solution in a tile t from 
the values given along the boundary of t. Let the solution of the boundary value 
problem 

(4.3) A^o = mo,i^t 

with boundary conditions 

{/a(xo), a;o G dta 
fbixo), Xo G dtb , 
fc{xo), Xo G dtc 

where the three sides of tile t are 9ta, dtb and dtc- Then, for ^ e i, we get 
(4.5) 

$(C) = j Qt{i,xo)fa{xo)dxo+ j Qt{ixo)fb{xo)dxo+ J Gt{£,,xo)fc{xo)dxo,£, et. 

dta atb ate 

Since V consists of copies of t, the solution of Eq. ()3.10p in V is the sum of integrals 
like Eq. ()4.5p . The A eigenvalue of operator A is called degenerate if in the case 
fa{£,) = fbiO = /c(C) = there exists a not identically zero <i>(f) solution over a 
tile t. The degenerate eigenfunctions are related to individual solutions over a tile 
with zero values fixed on the boundary. 

Corollary 4.1. For non- degenerate eigenvalues X, there is a one-to-one map be- 
tween the solution <I'(^),^ G t and the conditions prescribed on the boundary dt. 
When the conditions ifai{xo), fbiixo), fci{xo)) and {f a2{xo), fb2{xo), fc2{xo)) given 
along the respective boundaries of ti and t2 are linearly independent, the corre- 
sponding solutions ^i{S.),^2i£,) are also linearly independent. 

Corollary 14.11 connects the solution in a DiV V to the solutions along the bound- 
aries in V. When dealing with a Cauchy boundary value problem, we have only to 
deal with the internal boundaries. Let us assume that there are K internal sides (c.f. 
Definition l2.ip in V and pretend the solution fk{£,) to be known along internal side 
k. We use the vector form of $(a;), in which N functions F(^) = ($i(0, . . ., $Ar(^)), 
are encountered, one from each copy of t. Then, 

(4.6) FiO = Qv{0, 
where elements of w(^) are: 

(4.7) VkiO = J gti^,xo)fk{xo)dxo,k = l,...,K. 

Otfc 

Q is an i^T X matrisH, and assigns the K internal sides to the N copies of t. 
As we see, expression (|4.6p has two components, f (x) depends only on tile t and 
operator A, hence we call it the physical part of the solution. From now on, v,Vi, 
V2, . . ■ denote K tuples, underlined letters denote A^-tuples. Because of (|4.ip . and 
the assumed linearity of A, we get 

(4.8) AFiO = AQviO = QAviO = AF(e) 



Since every reflection creates an internal side, K=N-1. 
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for any ^ E t. Note that v depends also on the considered eigenvalue. The formal 
solution ()4.6|) reports us that the source of any non-idcntically zero solution derives 
from at least one non-identically zero function given along some boundary. We are 
dealing with a homogeneous problem therefore the solution is taken as zero at the 
external boundary of V. When the solution is not identically zero, we may assume 
that it differs from zero along the internal boundaries. In general, each boundary 
may be independent and the solution to boundary value problem (|4.3p must be 
linear combination of the solutions coming from the internal boundaries. Since the 
Green's function preserves linear dependence and independence, the dimension of 
the solution will not exceed the number of internal sides. 

On the other hand, Q depends only on the structure of V, hence we call it 
the structural part of the solution. Pattern (j4.6|) can be used to derive an integral 
equation set for the solution along the internal boundaries and may serve as basis 
for approximate solution methods. Now we are interested in the structural part. 

Lemma 4.2. The auxiliary matrix Xof volume V and the structural matrix Q in 
Eq. ( [4.6p are related as X~QQ^ . 

Proof: By the proposition, matrix element Xij is the dot product of the i*'' and 
j*^ rows of matrix Q. Hence, Xu equals the number of internal sides of copy i. 
Element Xij is 1 if copy i and j share an internal side, and zero otherwise. This is 
just the definition of X. □ 

Now we address the question of two discretized volumes Vi and V2 being equis- 
pectral. 

Theorem 4.3. Let discretized volumes volumes Vi and V2 such that 

(1) Vi and V2 are composed of the same tile t. 

(2) In Vi and V2 the number of copies (N) of t are equal. 

(3) Along the external boundary of Vi and V2 the number of sides are equal by 
side types. 

(4) It follows from item 3 that along the internal boundaries of Vi and V2 the 
number of sides are equal by side types. 

Let the formal solutions over Vi and V2 be 

Let the component v of the formal solution such that 

AVjn = \Vra,m =1,2 

holds. Furthermore, let the auxiliary matrices Xi and X2 be similar. Then, and 
only then Vi and V2 are equispectral. 

Proof: With appropriate norm, the eigenvalues are given by 

(4.9) A,„ = {AQ„^v^{0; Q^Vr,m),m = 1,2. 

Under the stipulated conditions Qi, Q2 are NxK matrices, where N is the number 
of copies of tile t in Vm , and K is the number of internal sides in either Vm ■ 

First we show that if the auxiliary matrices Xi and X2 are similar, then Vi and 
V2 are equispectral. Using (|4.9p . we get 

Ai - A2 = (Aui; QtQiVi) - {AV2; QtQ2'"2)- 
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Because of the similarity of the auxiliary matrices, we may write Q2 U Qi Qi U'^ , 
and V2 = Uvi where f/ is a unitary K x K matrix. Now we perform a series of 
elementary, identical, and algebraic transformations: 

(4.10) (Avr, Qt Q2V2) = (AUvi; Q+ Q^Uvi) = 

(4.11) ( UAvi; Q+ Q2 Uvi) ^ ( UAvi; UQt Qi U'^ Uvi) = 

(4.12) = {Avi;Q+Q^Vi). 

In (|4.10p . we expressed V2 by Vi, in (|4.1ip we expressed Q2Q2 by Qx and 
exploited that operator A acts only on Vi, and in ()4.12p we used the properties of 
the dot product. Using the last form, we get 

Ai - A2 - {Avx-{QtQxVi) - QtQiVi)^Q, 

from which follows Ai = A2. 

Now we show that if Ai — A2 then the auxiliary matrices are similar. By the 
assumption, we have 

[Avi] QtQiVi) = iAv2; QtQ2V2)- 

From that expression, operator A may be left out as it brings in a multiplication 
by the same number. Thus we have 

(vi; QtQlVl) = {V2] QtQ2'"2)- 

That expression holds for any components of vectors Vi and V2 because in the 
general case the solution along the K internal faces are independent, and the Green's 
function preserve that. Then, the two components of the dot product may differ 
only in a rotation: V2 — Uvi and Q2 Q2'"'2 — UQ^ Q^vi. From this immediately 
follows Q2 Q2 — UQi Qi iJ~^. We have to remember that Q2 Q2 is not the matrix 
X2, and from the similarity of Q2 Q2 and Qf Qi does not follow immediately the 
similarity of matrices X2 and Xi. We need some matrix theory |26j to go on. 
The N X K structural matrix Qi can be written as 

(4.13) UDV^ 

where U is an N x N orthogonal matrix, D is N x K matrix, its first K rows form 
a diagonal matrix, and the K + 1, . . . , N rows contain only zeros, V is an K x K 
orthogonal matrix. Since 

(4.14) QQ+ = UDV+VDU^ ^ UDDU^ , 

it follows that the nonzero eigenvalues oi QQ^ and Q are the same, therefore 
Xi and X2 are similar matrices as stated. □ 

Let us return to the case of a degenerate eigenvalue. When A is a degener- 
ate eigenvalue, we get a non-zero solution on a tile with identically zero values 
prescribed along the boundary of t. With a triangular tile t, the copies of t in 
discretized volume V fall into two categories (say black and white) as we can color 
V by the black and white colors. Let u; be an TV tuple with elements -1-1 or -1 in 
position i when copy ti is black or white, respectively. 

Lemma 4.4. The auxiliary matrix of the above considered discretized volume V, 
has the following property: Xw=0. The solution F{^) of eigenvalue problem I13.10\) 
has the following property: either F{^)=^{^)w where $(^) is the solution of problem 
fgJO)) on the tile t, or F{^)w=0. 
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Proof: The first part of the statement is an immediate consequence of the struc- 
ture of the auxihary matrix. The second part of the statement is a particular case 
of Hersch's theorem fTT . □ 

As elements in w are of opposite sign on adjacent tiles, F_{£,) = ^{0 * ML implies 
values for ^ on the inner boundaries. If v, (the physical part of the solution) 
does not vanish on internal sides, there is a vanishing linear combination {F_{S,)w) 
of solution values on respective points of different tiles. 

Theorem 4.5. Let F_„i{£,) = QmVm{0 ^ solution to Eq. iS.lO]} over DiV Vm, 
for ^ ei, m—1,2 with appropriate boundary conditions at the boundary ofVm- Let 
^mWm=0 for m=l,2. If there exists a matrix M such that MFi{S.)—F2{S,), then 
maps vector W2 wi ■ 

Proof: Since F_^{^) is a nontrivial solution, Vm{C} is not identically zero, in 
LemmaHm F_m{C} = '^{C)MLm is excluded, for m = 1, 2, we have u;^-F^(^) —Q^m — 
1,2. Since MF^{i) = £2(0, and w^£i(0 = 0, we have MwtF^{i) = 0. On the 
other hand, ^^£2(0 = ^2^-^£!i(0 = = ^i^Z!i(C)i from which immediately 
follows w_i = M^w_2- Q 

Theorem 14.51 is a strong constraint on the possible matrices transforming the 
solution of one boundary value problem into the solution of the other boundary 
value problem. 

4.2. Transplantation rule. Assume that discretized volumes Vx and V2 are eq- 
uispectral. What is the connection between the solutions $i(a;) and $2(2;) or, in 
vector form, between £i(^) and Fj2{^)l In accordance with Theorem 14. 3[ we may 
put the eigenfunctions into vector forms: ^^(a^) = £m(f )i "rn = \,2. Furthermore, 

(4.15) £„(0= g™^^m(0,"^ = l,2; 
and here 

(4.16) V2 = Uvi. 

This permits one to write down the matrix transforming the eigenvectors into each 
other. 

Lemma 4.6 (Transplantation of the solutions). The vector forms of the solutions 
are connected by the following linear transformation: 

(4.17) £2(0 = Q2U[QtQi]-'QtF,iO - MF.iO. 

Proof: From (|4.15p we get 

(4.18) QtF^iO 
and 

(4.19) QtF,iO^[QtQ^]v^iO 

Now the determination of M is straightforward: 

F^iO - Q2V2 - Q2UV, - Q2U{Q+Q,y' Q+F, 

where the first equality is the definition of £2, there we utilized (|4.16p . and expressed 
vi from (j4l9l) . □ 



12 



MIKLOS ANTAL AND MIHALY MAKAI 



Lemma 4.7. The transplantation matrix M is the solution of the following equa- 
tion: 

(4.20) M= Q2[QtQ2r'QtMQ,[QtQ,]-'Qt 
Proof: From (I4l8ll we get V2 = [Q^Qa]" QtE.2, 

using the transplantation 

matrix M, we get V2 = [Q2 Q2] ^ Q2 -^^^i- Here we use the definition of F_j^: 
^2 = [Q2 Q2] ^ Q2 MQiVi. We obtained an expression connecting Vi and V2, this 
is just matrix U: U= [Q2 Q2] ^ Qt ^Qi- Now we substitute this expression into 
the second equation in (|4.17p and arrive at (|4.20p . □ 

Lemma 4.8. When the transplantation matrix M is known, the transformation of 
the solutions along the internal sides is given by 

(4.21) U=[Q+Q2]-'QtMQ,. 

Proof: When F 2 = MF_^ and V2 Uvi, using (|4T8ll we get ^2 = [Qt Q2] Q2E-: 
In the last term we substitute F_2 ~ MF_i = MQ^Vi and the claim follows. □ 

4.3. Constructing equispectral volumes. The present Section is devoted to 
the problem of finding equispectral volumes. Our analysis is based on G and 
r, the group and graph associated with discretized volume V. Let us start by 
the investigations initiated by Sunada [T3]-[13. Let gSG and let {g} denote the 
conjugacy class of g in G. Let Gi, G2 CG subgroups in G. We say (G, Gi, G2) 
to form a Sunada triple if the number of elements from subgroups Gi and G2 are 
the same in every {g}. Let M be a manifold. 

Theorem 4.9 (Sunada theorem). Let V be a compact Riemannian manifold, G 
a finite group acting on M by isometrics. Suppose that (G, Gi, G2) is a Sunada 
triple, and that Gi and G2 act freely on M. Then the quotient manifolds Mi=Gi 
\M and M2 =G2 \M are equispectral. 

Starting out from the Sunada theorem, Gordon and Webb showed how to 
construct planar regions in which the Laplace operator is equispectral. 

One can find discretized volumes Vm equispectral to Vq by the Sunada theorem 
so that one searches for Sunada triples. To this end, we have to find in the group Gq 
associated with Vq subgroups of order N. The GAP program [21] offers means to 
solve that task. In Appendix A, we show such an algorithrrQ. The result is a coset 
representation of Go in terms of subgroups Gi and G2, from that the construction 
of Vi=Vo and V2 is straightforward. 

Following the recipe by Gordon and Webb [TQ , we study the graphs associated 
with Vo and V^. In that, the following conjecture is utilized. 

Conjecture 4.10. Let Vi and V2 be equispectral volumes with associated graphs 
Fi and F2 and associated groups Gi and G2. Lf graphs Fi and F2 with respective 
edges Sj and Vi, are isomorphic, and the isomorphism of the edges is Si — giVig^^ , 
then Vi ~ V2 does not hold unless all g^ are automorphisms of tile t when Gi and 
G2 are isomorphic groups. 

Buser |14| has presented equispectral graphs which are not isomorphic but the 
associated discretized volumes are not planar ones. Buser's graphs are derived from 



^ Courtesy of Dr. Erzsebet Lukacs of BME Institute of Mathematics. 
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abstract groups and the adjacency proposed by the coset representation exclude the 
associated geometry to be planar. Buser, Conway, and Doyle [17] presented planar 
equispectral graphs, their "warped propellers" are equivalent but if the underlying 
tile t is not a regular triangle, those are really equispectral. By means of the tools 
of projective geometry, several further equispectral pairs must be created. 

The equispectral DiVs seem to be related to isomorphic graphs and isomorphic 
groups. This is formulated in conjecture 14.101 In Appendix B, we give a proof for 
a special case of that conjecture. 



5. An example 

Since isomorphic graphs contain the same number of vertices of degrees 1, 2, 
and 3, one seeks equispectral discretized volumes by specific features (length of 
walks, number of vertices of degree 3) of the graph. To this end, we have analyzed 
V composed of seven regular triangles. There are 25 non-isomorphic graphs, the 
order of the associated groups are 5040=7! for 12 graphs, 2520 for 10 graphs, and 
168 for 3 graphs. In accordance with conjecture I4.10i equispectral volumes may 
exist within a family of non-isomorphic graphs provided the tile t is not symmetric. 
We show an example first presented in Ref. 13J. The order of either associated 
group is 2520, the generators of the group Gi assciated with the DiV on the left are 
ci = (1, 2) (5, 6), 6; — (3, 5) (2, 4), and a; = (4, 6) (5, 7), whereas the generators of the 
group Gr associated with the DiV on the right are C2 = (1, 2) (5, 6), 6^ — (3, 5) (2, 4), 
ttr = (3, 7)(2, 6). The generators are related as Cr — q, br = bi, and — aibiai. The 
associated graphs are l-{c}-2-{b}-4-{a}-6-{c}-5(-{b}-3)-{a}-7 (left) and 7-{a'}-3- 
{b'}-5-{c'}-6-{a'}-2-(-{c'}-l)-{b'}-4 (on the right). (Here i-{c}-j stands for vertices 
(i,j) connected by an edge of type c. A vertex of degree 3 has two connecting walks, 
one of them is put into parentheses). The discretized volumes shown in Fig. 4, are 
the so called Schreier graphs of groups Gi and G2. These graphs are associated 
with the discretized volumes depicted in[3l The DiV on the left/right is associated 
with subscript l/r, respectively. The auxiliary matrices are 



(5.1) 



X, = 



/ 1 1 \ 

1 2 1 

10 10 

1 2 1 

1 3 1 1 

1 1 2 

V 1 1 / 



and 



(5.2) 



/ 1 1 \ 

10 10 

2 10 1 

10 10 

1 2 1 

10 12 

V 1 1 y 



1 3 
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rank{Xi) = rank{Xr) = 6, the eigenvalues are (0, 0.225377, 1, 1, 2.18589, 3.36041, 4.22833). 
The structure matrices are 



(5.3) 



Qi 



/ 1 \ 

1 1 1 

1 1 

110 
110 

V 1 y 



and 



(5.4) 



Qr 



/ 1 \ 

1 1 

1 
1 1 
1110 
110 

\ 1 y 



The first trial is to express the solution along internal boundaries as linear combi- 
nations of the solution along the corresponding internal boundaries. Since solutions 
along the side of the same type may be combined only, we may trj0 



(5.5) 



U 



( 



1 
10 
10 10 
1 
-10 10 
10 



-1 \ 




1 


0/ 



Using Eq. (|4.17p . we get 
(5.6) 



M- 



I 1 

1 



-1 1 

1 

1 

V 



\ 


-1 



1 

-1 

/ 



V 



1 

1 

-1 1 



1 
1 

-1 1 



-1 \ 

1 



-1 




1 / 



The two linearly independent matrix gives two variants of the transplantation ma- 
trix M.. We investigated if it is by chance. From the condition that the solutions 



^For simplicity's sake the normalization factor is omitted. 
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should be in C^, after a long computation one obtains the following general struc- 
ture of transplantation matrix M: 

I 



(5.7) 



M 



a 


b 


a 


^b 


—a 


-b 


—a 


\ 


b 


a 


—a 


b 


-b 


—a 


a 




a 


a 


b 


b 


a 


—a 


-6 




b 


b 


-b 


a 


—a 


a 


—a 




a 


b 


a 


a 


b 


a 


b 




b 


a 


—a 


—a 


a 


b 


-b 




a 


a 


b 


—a 


-b 


b 


a 


) 



V - 

Here a and b are independent, using the choices a — 0,5 = 1 and a ~ l,b = 
one gets the two linearly independent transplantation matrices. Gordon, Makover 
and Webb [25] state that the above decomposition corresponds to decomposing 
a representation of the group G/=Gr by 7 x 7 matrices into two, nonequivalent 
irreducible components. The degenerate solutions are related to 

(5.8) 
(5.9) 



- (1,-1,-1,1,1,-1,-1) 
w. = (1,-1,1,1,-1,1,-1). 



In the analysis of discretized volumes of practical calculations, the presented 
considerations may serve only as a theoretical investigation for two reasons. The 
first one is the large number of copies in a discretized volume. The algorithm 
presented in Appendix A works only for N<10, which is really far from the practical 
applications (N>1000). The second one is in the limited choice of equispectral 
volumes: they have the same number of tiles, the external and internal boundary 
types should correspond in the two DiVs. Some of the limitations may be mitigated 
or eliminated in the future. 



6. Concluding remarks 

The subject of the present work is the discretized volumes (DiVs), a tool fre- 
quently encountered in applied mathematics. Although our investigations have 
been limited to DiVs composed of a small number of copies of a tile, the conclu- 
sions are applicable to real problems as well. We considered the solution of an 
eigenvalue problem in a DiV, in which the involved operator is linear, its symmetry 
include the symmetry group of the plane M^, and the solution belongs to C^. The 
goals of the present work are: 

• To find algebraic descriptions for discretized volumes (DiVs) in order to 
find out whether two DiVs are equispectral or not. 

• To find criteria for two DiVs to be equispectral. 

• To find a transplantation recipe with which we are able to transform the 
solution from one DiV into another one. 

Lemma 4.2 connects a component of the formal solution, the structural matrix Q 
and a description of the DiV, the auxiliary matrix X. Theorem 4.3 gives sufficient 
and necessary conditions for two DiVs to be equispectral. The theorem uses only the 
structural part of the formal solution, therefore the theorem applies to a large class 
of operators involved in the eigenvalue problem. The eigenvalue connected with zero 
solution along every internal sides is a degenerate case: that degenerate problem 
should be discussed separately. This is done in Theorem 4.5. The transplantation 
rule has two components: the transformation of the solution along internal sides 
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and the transformation of the solution along the copies if the tile. The relationship 
among them is established by lemmas 4.6-4.8. As to the graphs associated with 
a DiV, we formulated conjecture 4.10 in the hope to attract the attention of the 
experts, perhaps algebraic geometry means prove more successful here. 
We achieved the following results: 

(1) Formulated an eigenvalue problem with Cauchy boundary condition for 
DiVs and associated algebraic descriptions with the DiVs, viz. a graph and 
a group. 

(2) By analyzing that group and graph, one can establish relations between 
eigenvalue problems on specific DiVs. 

(3) We improved the formal solution given in Ref . [22j . The improvements have 
lead to Lemma l4^ and Theorems l4.3l and l4.5l We formulated necessary and 
sufficient conditions for two DiVs to be equispectral: when the eigenvalues 
of the respective auxiliary matrices are the same then and only then are 
equispectral two DiVs. 

(4) We pointed out a surprising relationship between the eigenvalues of the aux- 
iliary matrix Xand the spectrum of operator A over DiV V. Furthermore, 
the degenerate eigenvalues of V are associated with the zero eigenvalue of 
X. Equispectral DiVs Vi and V2 have similar auxiliary matrices. 

(5) Conjecture 14.101 leads to a classification of possible DiVs. We made the 
classification of the DiVs composed of seven triangles, the most frequently 
encountered problem in the literature [13]- [19], [22]. Known equispectral 
DiVs confirm our observations. 

On the other hand, our investigation addresses new problems as well. 

• It is a question if Cayley graphs of planar equispectral volumes are always 
isomorphic (Conjecture 14. 10|) . 

• Whether there is an algorithm to solve the problem considerably faster than 
with the algorithm in Appendix A? Allows that speed for solving practical 
problems (N>1000)? 

• Buser. Conway, Doyle, and Semmler expressed their guess: equispectral 
volumes are rather scarce. Is it true for very large DiVs (i.e. for NwlO^)? 

• Can the formal solution (|4.6|) be generalized for DiVs with different tiles? 

Discretized volumes (DiVs) offer a way in which non-equivalent geometries can 
be found to solve an equation and a simple transplantation rule can also be given. 
Then the transplantation is exact. It is true that in the present form the procedure is 
simple and of restricted use. However, there are reserves to be exploited. Triangles 
that can be transformed into each other by a linear map may allow for an extension 
of the method. 

7. Appendix A 

The GAP program given below has been written to find discretized volumes Vi 
equispectral to a given discretized volume Vq. Vq is defined by its three generators. 
The discretized volumes Vi are characterized by the number of copies connected by 
sides of type a, b and c. The version presented here has been made to verify if the 
program finds the equivalent discretized volume pairs discovered by Gordon and 
Webb, see Fig. 3. 

# # Discretized volume V is given by the three generators of group # 
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G. # Now V is made up from 7 copies ot tile t # # The search starts 
from VO given also by three generators d, e and f 

\# d:=(3,7)(2,6); e:=(2,4) (3,5) ; f : =(1 , 2) (5 , 6) ; {\#} LogTo("Tri7") ; 
G:=Group(d,e,f ) ; Size(G); CT:=[];; ri:=[]; 
cc : =ConjugacyClassesSubgroups (G) ; ; 
ccr:=List(cc,x-$>$Representative(x)) ; ; 

ccf :=Filtered(ccr,x-$>$Size(x)*7=Size(G)) ; ; {\#} {\#} In ccf, we 
have collected all subgroups of index 7 {\#} Next we apply 
generators d, e, f, to them {\#} for i in [1 .. Size (ccf ) ] do 
Append (CT, CosetTableBySubgroup (G,ccf [i] )) ;Append(ri, [i] ) ; od; {\#} 
{\#} In CT, we have the discretized volumes represented by coset 
tables 

{\#} (One line is associated with ech generator and its 

inverse) 

{\#} 

8. Appendix B: A proof for a special case of Conjecture 14.101 

Proposition 8.1. A Cayley-graph is given (T^), which is depicted with isosceles 
triangle nodes (N), the graphic is denoted by A. If the equal sides of the node are 
a and c, and the third side is denoted by h, then a change of a and c in the graph 
means a reflection of the graphic: TA = B. (The graph we get after the change is 
r^. The reflection is denoted by T, and B is the graphic we get after the change 
in the graph.) 

Proof: For the proof we use complete induction. If we start drawing and 
at the same node (A'o) fixed at a position, then the statement holds for this 
zeroth order case: the node is symmetric to the altitude perpendicular to the base: 
TMp = Mq. Side a and c are mirror images of each other, side b is reflected to 
itself. 

Let us denote the sides of triangle Ni by ai,bi,Ci. Thus: Toq — Cq .Th^ = 

Then, in each step a new triangle (Ni) is obtained and drawn in both graphics by 
reflecting a formerly got triangle. We show, that if the statement is true for step n 
(that means TA„ = i?„), then it holds for step {n + 1) (meaning TAn+i — 
We distinguish two cases: 

If reflection r (here r denotes a reflection to a side and it affects only one triangle, 
when T is the reflection to the ordinary axis of symmetry and it affects the whole 
graphic) is made with respect to side fo, then we are reflecting to such sides in the 
two graphics, that are mirror images of each other: Tbf ^bf . It also holds for the 
actually reflected triangles: TNf = Nf (as TAi = Bi Vi < n-re). Using Lemma 
18.21 the statement concludes. It is subsequently true, that: Taf = cf,Tbf — 

If the reflection is made with respect to side a or c, then beside TA^^ = Nf we 
can also observe that Ta^ = Ci and Ta = Oj, so Lemma \K2[ can be used as well. 
So: Taf = cf , Tbf = bf , Tcf = af Vi. 

Lemma 8.2. Let us consider a reflection with respect to a certain (a, b or c) side 
(t). IfTNf = Af and Tr^ = r^ then T{rN^) = rNf , where r^ and r^ denote 
the reflection axis in graphic A and B respectively. 
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Proof: the reflection of a triangle can be described with the reflection of its 
vertices. Two vertices of are on the reflection axis, so they are transformed by T 
into the appropriate vertices in the other graphic. Third vertices are mirror images 
of each other, as it can be seen from basic coordinate geometrical considerations. 
(Reflection r should be considered as a coordinate geometrical transformation, then 
it is trivial, that T transforms those vertices into each other.) So: Taf = cf , Tbf = 
bf,Tcf=af. □ 
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Figure 1. A discretized volume from four triangles 
The sides are labelled by the generators 




Figure 2. Sequence of gluing 
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